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Abstract 

Taking into account the structural transition and long-range interaction (lattice effect), we resort to the 
Ewald-Komfeld formulation and developed Maxwell-Gamett theory for uniaxially anisotropic suspensions 
to calculate the effective permeability of inverse ferrofluids. And we also consider the effect of volume 
fraction to the magnetophoretic force on the nonmagnetic spherical particles submerged in ferrofluids in 
the presence of nonuniform magnetic field. We find that the coupling of ac and dc field case can lead to 
fundamental and third harmonic response in the effective magnetophoresis and changing the aspect ratio in 
both prolate and oblate particles can alter the harmonic and nonharmonic response and cause the magne- 
tophoretic force vanish. 
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I. INTRODUCTION 



Magnetic particles are used to label and manipulate biomaterials such as cells, enzymes, pro- 
teins, DNA, to transport therapeutic drugs in the applications such as bioseparation, immunoas- 
says, drug delivery, and to separate red and white blood cells yj]. Magnetizable particles expe- 
rience a force in a nonuniform magnetic field caused by the magnetic polarization. This phe- 
nomenon, called magnetophoresis (MAP) has been exploited in a variety of industrial and com- 
mercial process for separation and beneficiation of solids suspended in liquids [i2fl. Virtually, 
almost all practical magnetic separation technologies exploit the magnetophoretic forces. Such as, 
belt and drum devices, high-gradient magnetic separation (HGMS) systems [3] etc.. 

In recent years, inverse ferrofluids with nonmagnetic colloidal particles (or magnetic holes) 
suspended in ferrofluids have drawn much attention for their potential industrial and biomedical 
applications si 6[ 0]. The size of the nonmagnetic particle (such as polystyrene particle) is 
about l~100/im. Ferrofluids, also know as magnetic fluid, is a colloidal suspension of single- 
domain magnetic particles, with typical dimension around lOnm, dispersed in a liquid carrier 
Because the external magnetic field will induce magnetic dipole moment in the nonmagnetic 
particles, when the field intensity increases, the inverse ferrofluids will solidify to crystals due to 
dipole interactions, which is similar to ER and MR behavior. Since the nonmagnetic particles are 
much larger than the magnetic fluid nanoparticles, the magnetic fluid is treated as a one-component 
continuum with respect to the larger nonmagnetic particles [8,1^. Recently, Jian, Huang and Tao 
first theoretically suggested that the ground state of inverse ferrofluids is body-centered tetragonal 
(bet) latticejiol] , just as the ground states of electrorheological fluids and magnetorheological 



fluids lU 



1211 ■ Motivated by this study, we propose an alternative structural transition from the 



bet to the face-centered cubic (fee), through the application of external magnetic field. In this 
work, we will investigate the magnetophoretic force on the particles during the phase transition in 
the crystal system considering the effects mentioned above. 

Finite-frequency responses of nonlinear dielectric composite materials have attracted great at- 
tention in both research and industrial applications during the last two decades [|23|] . In particular, 
when a composite containing linear permeability particles embedded in a nonlinear permeability 
host medium is subjected to a sinusoidal alternating magnetic, the magnetic response w ill g ener- 



ally consist of frequencies of higher order harmonics under ac field 
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sob. The 



main aim of the present paper is to study the the effects of the geometric shape and volume fraction 
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of the particles on the nonlinear ac responses (harmonics) of the crystal system, then the magne- 
tophoresis of the particles, considering the local lattice effect. We shall use the Edward- Komf eld 



formulation 114 , 



15i 



1811 to derive the local magnetic fields and induced dipole moments in inverse 



ferrofluids and then perform the perturbation expansion method to obtain the fundamental and 
higher harmonics in two cases: single ac and dc-ac magnetic field. 

Then, we apply a nonuniform magnetic field to investigate the magnetophoresis of the nonmag- 
netic particles submerged in the ferrofluids, by taking into account the effect of structural transition 
and long-range interaction (lattice effect). And we also consider the effect of volume fraction to 
the magnetophoretic force exerted on the nonmagnetic particles. 

This paper is organized as follows. In Sec. II, we introduce the nonlinear characteristics in 
the ferrofluids and the theory of magnetophoresis. By using the Ewald-Komfeld formulation and 
the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, we cal- 
culate the effective permeability of the inverse ferrofluids for different harmonic response under 
time-varying magnetic field in Sec. III. Then we investigate the nonlinear response of magne- 
tophoretic force on the nonmagnetic particles in ferrofluids. This paper ends with a discussion and 
conclusion in Sec. IV. 
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II. THEORY AND FORMALISM 



A. Nonlinear characteristics in the ferrofluids 

In the standard literature, we often consider the magnetic induction B to be proportional with 
the field strength Hq, and their relation is written as B = /xHo, where /i is the linear permeability. 
However, the simple relation is only valid at weak intensities. When in the real situation, especially 
under a stronger field, nonlinearities are introduced as B = yuHo + ,^|//opHo + r]\Ho\^Ho + 
where ^,r] are the third-order and fifth-order nonlinear coefficient and 771-^01"^ ^ ^l-^^oP ^ A*- 
Here we assume the nonlinearity is not strong and consider only the lowest-order nonlinearity for 
simplicity. Thus after dropping the symbol of absolute value in the equation above, the nonlinear 
permeability jl is obtained: 

/i = /ii + eii^o (1) 

For inverse ferrofluids, the peameability fx of host medium ferrofluid can be derived in the form 
as /i = yUi + ^iHq when the magnetic particles in the ferrofluid is nonlinear to the external field, as 
analyzed below. As shown in Fig.(l), considering the magnetic particles with nonlinear jlp in the 
ferrofluid , the local magnetic field Hi inside the particle induced by external field Hq is give by 

Hi = Ho (2) 

Here fi j represents the host medium permeability in the ferrofluid, and the nonlinear perme- 
ability fip of the magnetic particle can be expressed as jlp = Hp + ^ Hf with ^ being the nonlinear 
coefficient. For the spherical shape particles, the effective permeability yU/ e of the whole ferrofluid 



can be determined by the Maxwell- Gamett theory [|22h 



where / denotes the volume fraction of the magnetic particles in ferrofluid. Applying Eq.(2) 
and (3) and performing the Taylor expansion by taking ^Hf as a small perturbation, we can obtain 



^^''^ (l-/K + (2 + /W ^ [{l-f)^^p + i2 + f)^,f]{^^p + 2^^fr^ 
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Hence the nonlinear permeability of ferrofluid can be expressed as 

/i = yU/,e = /ii + ^iHq (5) 

in which the constants fx, = A'/[(i+2/)/^p-2Cf-i)M/] ^ _ J,'Jn\( 4-9 ^2 + 0[e] We 

emphasize that when the magnetic particles form chains under strong external field, Eq. (3) is sup- 
posed to modify by introducing demagnetization factor |3^, but the expression of jl will remain 
unchanged by appropriately choosing fii and ,^1. 

B. Local Magnetic Field and Induced Dipole Moment 

When the inverse ferrofluids system forms crystal under external field, there are two dominant 
factors effecting its effective permeability: the local field of crystal structure and the geometrical 
shape of the particles which makes up of the crystal. Below we will first investigate the local field 
factor by performing an Ewald-Kornfeld formulation so that the structural transition and long- 
range interaction can be taken into account explicitly. Considering the ground state of inverse 
ferrofluids as a body-centered-tetragonal (bet) lattice, which can be regarded as a tetragonal lattice, 
plus a basis of two particles each of which is fixed with an induced point magnetic dipole at its 
center. One of the two particles is located at a corner and the other one at the body center of the 
tetragonal unit cell, as shown in Fig.(l). Its lattice constants are denoted by a{= h) = ^q^^ and 
c= along x{y) and z axes. In this case we identified a different structure transformation from 
the bet to fee to bee by changing the uniaxial lattice constant c under hard-sphere constraint. As q 
varies the volume of unit cell remains f c = • In this way, the degree of anisotropy of tetragonal 
lattice is measured by how q is deviated from unity and the uniaxial anisotropic axis is along the z 
axis. In particular, q = 0.87358, 1 and 23 represent the bet, bee and fee lattice, respectively. 

In case of an external magnetic field Hq along the z axis, the induced dipole moment m are 
perpendicular to the uniaxial anisotropic axis. The local field Hr at the lattice point R=0 can be 
determined by using the Ewald-Kornfeld formulation 



In this equation, 71 and 72 are two coefficients, given by, 

erfc(r/r) 2r] / 2 2n rn\ 
liyr) = o — + exp{-r] r ) (7) 
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3erftM + ,V ^ % ) ( V) (8) 

where erfc(r7r) is the complementary error function, and r/ is an adjustable parameter making the 
summations in real and reciprocal lattices converge rapidly. In Eq. (7), R and G denote the lattice 
vector and reciprocal lattice vector, respectively, 

R = C,{q~^lx + q~^my + qnz) (9) 



G = —{q^ux + q^vy + q ^wz) (10) 
where, Z, m, n, m, f , w are integers. In addition, Zj and Rj are respectively given by. 



j - 1 j - I 
Zj = I — , Rj = \R — {ax + ay + cz)\ (11) 



and n(G) = 1 + exp[i{u + v + w)/tt] is the structure factor. The local field will be computed 
by summing over all integer indices (j, /, m, n) ^ (1, 0, 0, 0) for the summation in real lattice 
and {u, V, w) ^ (0, 0, 0) for the summation in the reciprocal lattice. Here we consider an infinite 
lattice. For finite lattices, one must be careful about the effects of different boundary conditions 
I19I1 . Because of the exponential factors, we may impose an upper limit to the indices, i.e., all 
induces ranging from —L to L, where L is a positive integer. 
Now, let us define a local field factor a, 

a = (12) 

An m 

It is worth noticing that a is a function of single variable q. Also, there is a sum rule 2a± +a\\ 



=3 [2^|21I], where «|| denotes the local field factor in longitudinal field case, while a± denotes 
the local field in transverse case. Here the longitudinal (or transverse) field case corresponds to 
the fact the Hq field is parallel (or perpendicular) to the uniaxial anisotropic z axis. 



C. Nonlinear AC response in the Cases of Spheroids 

Now we are in the position to take into the account the geometrical shape effect of the particles 
by investigating nonspheres: prolate and oblate. In MR(ER or inverse ferrofluids) fluid, the shear 
stress and permeability(permittivity) will be enhanced by changing the ratio of the ellipsoidal 
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particles in the chain CSll such as p-phenylene-2, 6- benzobisthiazole polyparticle. This induced 
polarization can be described in demagnetizing factor, which includes the three principal axes of 
the spheroid. For a prolate spheroid(a = b < c), the longitudinal demagnetizing factor gi along 
the major axis c is given by gi 



fcln(k+v^k^) 

(fc2-i)i 



fc2-l 



with A; = ^; for oblate case(a < b = c). 
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fc2 arcsin V^'^^"^ 

1 ^f^^ — ^ 1^ ■ There is a sum rule for the longitude and transverse demagnetizing 



factor [122, 



3211 that gi + 2gt = 1. For prolate spheroid, because it can easily align with its major 



axis(c) along the direction of the external field, while the dipole moments of nonmagnetic particles 
are oriented opposing the external field, here we only focus on the longitude case. It is also worth 
noticing that the spherical shape of particles is included when gi = gt = \. 



N 



1122L 



ow We can use the developed Maxwell-Gamett theory for uniaxially anisotropic suspensions 



3l|] considering both the local lattice and geometrical shape effect to derive the effective 



permeability [l32ll 



Similar to the technique applying in Eq.(4), combining Eq.(l) we admit 



gi{fio - fi) 



(13) 



HilaHi + gijxa - 3p)y] 
+ 9ixy] 



afif + gfxy^[xa — 3p] + 2agifiixy 3(?f/iQ/i^pa;^^if, 



a[y9i - i^iVii^i + gixvY 



(14) 



with X = p — l,y = fii — fiQ. 



In the presence of external oscillating time- varying magnetic field cite jlTll . the magnetic 
particles will have nonlinear characteristics and the inverse ferrofluids system will induce har- 
monic response at different frequencies without phase transition lllin . In the experiment, the 
second term or higher order terms on the right-hand of Eq. (14) can be obtained using 
mixed-frequency measurements jisl . In this paper, we will focus on the characteristics of the har- 
monics of effective permeability for the inverse ferrofluid and their nature on the magnetophoresis 
when particles put in: i^o-, Hu:-, l^2w-- In case of an external sinusoidal magnetic field(ac) of the 
form 



-o[e 



Ho = HacSinuot (15) 

Applying Eq.(15) into Eq.(14), we can obtain zero, second and fourth order harmonics for the 
effective permeability: /Xe = /^o + /X2tjCos2a;t + /i4^cos4a;t. Here we give the expression for n^^: 
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r2 Sgffi^px 
'8a[gixy + giy] 

On the other hand, under the application of a dc and ac field, namely. 



Hq = Hdc + Hacsmut (17) 

Using the same technique we can rewrite the effective permeability as /ig = /io + /io^sincut + 
yU2a;Cos2ci;t + yUa^sinStut + yU4^cos4co't, for extracting different harmonics. 

D. Magnetophoresis in inverse ferrofluids system 

Now we will investigate the magnetophoretic behavior when particles are put in inverse fer- 
rofluids. In a standard case, a magnetically polarizable object will be trapped in a region of a fo- 
cused magnetic field, or nonmagnetic particles in a magnetizable liquid(ferro fluid), provided there 
is sufficient magnetic response to overcome thermal energy and the magnetophoretic force ||2|]. 
Here we consider the nonmagnetically homogenous particle in the nonlinear inverse ferrofluids 
for simplicity, and the magnetic field is assumed to be dc and hysteresis is ignored. For a magneti- 
cally linear particle under magnetophoresis, the effective magnetic dipole moment vector induced 
inside takes a form very similar to that for the effective dipole moment of dielectric paricle 161. 
m = ATiRlj^jbHo, where Rejf is the effective radius of spherical or spheroidal particle and b is the 
clausius-Mossotti factor(CMF) along the direction of external field. Here we focus on the magne- 
tophoresis of spherical particle in inverse ferrofluids, and CMP b is written as: b = ^'^^^^ where 
Ii2 is the permeability of nonmagnetic particle, and jie is the permeability of the host medium: 



inverse ferrofluids. It should be noted, that in a general case, for spheroidal particle, CMF b is 

1 A'2-A'e 
3 Me+gi(M2-Me) ■ 



written as \ J^'^, ^1 — The force exerted by a nonuniform magnetic field on the dipole can be 



obtained: 

— * 

Fdipole = fJ'em ■ VHq (18) 

Combining the effective magnetic dipole moment, the magnetophoretic force exerted on the 
spherical particle in a nonuniform magnetic field can be written as, 

Fmap = 2Tt^ieR'bVHl (19) 
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Thus the magnetophoretic force on particles are proportional to the particle volume, the polar- 
ization difference and distribution of the geometric field gradient. There are two kinds of magne- 
tophoresis depended upon the relative magnitudes of ^2 and /ig. Particles are attracted to magnetic 
field intensity maxima and repelled to the field-generating electrodes when 6 > as positive mag- 
netophoresis and negative magnetophoresis corresponds to 6 < 0. Based on the harmonics of 
effective permeability, the expressions for the CMF can explicitly be obtained when the lattice 
structure of the host medium is changed, as shown in Fig. 2. The degree of anisotropy q plays 
an important role in determining the local field factor a. There is a plateau during the structure 
transformation on the relationship between q and a [|29D, accordingly we can also predict similar 
behavior in CMF, especially in Fig. 2(b). It is found that increasing q causes both zero and fourth 
order harmonic CMF decreases in Fig. 2(a) and Fig. 2(c). Note that when q changes from 0.6 to 
0.8 all CMF have great decreasing and CMF in Fig. 2(c) changes its sign when the nonlinearity is 
the highest, thus demonstrating the sensitivity of CMF when the system is under structure changes. 
Although the fourth-order harmonics of the local magnetic field and induced dipole moment seems 
to be complicated, since the strength of the nonlinear polarization of nonlinear materials can be 
reflected in the magnitude of the high-order harmonics of local magnetic field and dipole mo- 
ment, the effective nonlinear part in the effective nonlinear magnetic constant can be expected to 
determine the magnitude of the fourth-order harmonics. 



ni. NUMERICAL RESULTS 



In the preceding section we will calculate the effect of spheroidal shape and volume fraction of 
nonmagnetic particles which forms bet lattice structure in the ferro fluids, on the harmonic response 
of CMF. For a special case, when we would like to conclude the nonharmonic response situation, 
the external field Hq is constant and not vary with time namely, it can be easily obtained without 
tedious calculation by setting cut = | for ac field in Eq. (15). Thus only the zero harmonic 
contributes. Similar conclusions can be obtained that only hQ + h^ — b^^ contributes in a dc and ac 
field case. 

In Fig. 3 and Fig. 4 we show the harmonic response of CMF for different shapes in ac field: 
prolate and oblate. The zero(or nonharmonic) response of CMF exhibits strong sensitivity to the 
prolate demagnetization factor(the shape) and some peaks are obtained, while for oblate CMF bo 
have slight change compared(Fig. 4). This can be briefly explained through self-consistent ap- 
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proach 113 311 by the spectral representation separating the material parameter s = {1 — hq/ fi) from 
the geometrical parameters. Without considering the crystal lattice effect, the effective permeabil- 
ity can be derived as /ie = — g-g;('i-p) ]- It is clear that decreasing gi leads to increasing /ig, 
while for prolate(oblate) case, increasing aspect ratio k causes increasing(decreasing) of gi. In 
the meantime, because of the existence of weak nonlinearity, the nonlinear polarization has just 
a perturbation effect, which can be neglected when compared to the linear part. Furthermore, we 
can predict that one or two shape-dominant crossover ratio at which there is no net force on the 
particle in the magnetophoresis. In Fig. 3(a), it is interesting that the crossover ratio is mono- 
tonically increasing as the nonlinearity increases. Fig. 4 shows that, as nonlinear characteristics 
increase, the harmonics response of CMF are caused to increase accordingly for oblate particles. 
We emphasize to point out that, although prolate and oblate particles can both exist in different 
system state, the alignment of one particle can only be along the longest axis(prolate) when the 
particle is homogenous under the external magnetic field. 

Fig. 5 and Fig. 6(a)-(e) show different harmonic response of CMF for prolate and oblate 
particle in dc and ac field. Comparing with in single ac field, the response nature of CMF keeps 
unchanged for the same harmonics. So the introduction of coupling between dc and ac field add 
more odd harmonic response. As discussed above the nonharmonic response of CMF can be 
obtained as + h^^ — b^^j, which is shown in Fig. 5(f) and Fig. 6(f). It is also found that lower 
harmonic response of CMF can be several orders of magnitude larger than the higher one. 

Fig. 7 displays for prolate, oblate spheroidal and spherical particle, the effect of volume fraction 
of magnetic particles in the ferrofluids on the harmonic response of CMF. It is observed from 
Fig. 7(a) that the nonlinearity for nonharmonic response has little effect on CMF. The reason is 
that in the system of interest only ferrofluids are assumed to be nonlinear while the particles are 
linear. As the volume fraction of particles increases, the nonlinear component within the system 
have slight change naturally, yielding unchanged CMF accordingly. The second and fourth order 
response of CMF in Fig. 7(b) and (c) shows great variety when the volume fraction changes. This 
enhancement of the harmonics may arise from the mutual interaction between the dipole moments: 
when increasing the volume fraction, the mutual interaction becomes strong and leads to shift of 
harmonic response. It is observed that for prolate particles, the nature of zero-order response of 
CMF is greatly different from that of oblate spheroidal and spherical particles, comparing Fig. 
7(a) with Fig. 7(d) and Fig. 7(g). 
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IV. DISCUSSION AND CONCLUSION 



In this work, we have investigated both the effect of lattice structure of crystal and the geomet- 
rical shape of particles embedded in the ferrofluids on the nonlinear response of magnetophoresis 
using the perturbation method. The nonlinear response of CMF under the structure transition is 
studied in detail. It is thus possible to perform a real-time monitoring of structure transforma- 
tion. We find that the coupling of ac and dc field case can lead to fundamental and third har- 
monic response in the effective permeability constant and the CMF respectively. We also find that 
change the aspect ratio in both prolate and oblate particles can alter the harmonic and nonhar- 
monic response and thereby causing the magnetophoretic force vanish. Furthermore, by taking 
into account the local-field effect from the mutual interaction, we examine the CMF response in 
various volume fractions, and find that increasing volume fraction can result in the enhancement of 
the strengths and the difference between the different nonlinearities and characteristic harmonics. 
Thus, it is possible to obtain a good agreement between theoretical predictions and experimental 
data by suitable adjustment of both the geometric parameters(for example, the particle shape) and 
the physical parameters(for example, the temperature, component ratio in the system, saturation 
magnetization, initial magnetic susceptibility). On the contrary, such fitting is useful to obtain the 
relevant physical information of magnetophoresis in inverse ferrofluid. 

We demonstate theoretically the shape effect of the effective magnetic constant in the inverse 
ferrofluids system. It is in fact that the magnetic particles in the ferrofluid are inhomogeneous 
in the real case. A first-principle study of the dielectric dispersion of inhomogeneous suspension 



is presented 1291 13411 . and can be used to analyze the magnetic dispersion in the dilute limit, in 



which the permeability of the ferrofluid is dependent on the frequencies of external magnetic field. 



Also when the inhomogeneity is simp 



ified as the shelled spheroidal particles in the ferrofluid 



system BSll . a self-consistent method Oal can be applied for the three-component composites. 



The magnetic anisotropy as well as geometric anisotropy can both be taken into account, and will 
be reported elsewhere. 

The nonlinearity under our consideration is weak, which is common in real situations under 
moderate fields, while in case of strong nonlinearity, the perturbation approach is no longer valid 
and should be adopted using self-consistent method instead [|37|l . And in the above discussion, we 
consider the longitude case under the external field, it can also be found constricted by the relation 
gi + 2gt = 1, the harmonic response of CMF have opposite behavior. 
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Figure Captions 

Fig.l. (Color online) Schematic graph showing the inverse ferrofluid system, in which the crys- 
tal consisted of nonmagnetic particles(permeability is denoted by hq) embed in the ferrofluid(/io). 
Ferrofluids contain nanosize ferromagnetic particles(/ip) dispersed in a carrier fluid(^/). 

Fig.2. (Color online) The zero, second, fourth harmonic response of CMF versus degree of 
anisotropy q for different intrinsic nonlinear characteristics in the prolate case. The bet, bcc and 
fee lattices which are respectively related to g = 0.87358, 1.0,25 are shown(dot lines). Parameters: 
A*o = A«2 = = 2,p = 0.1, A; = 1.2. 

Fig.3. (Color online) Cases of prolate spheroidal particles: ac field. Zero, second, fourth 
harmonic response of CMF versus k — c/a for different intrinsic nonlinear characteristics. Pa- 
rameters: /Uo = = 1, yUi = 2, p = 0.1,0; = 0.95351. 

Fig.4. (Color online) Cases of oblate spheroidal particles: ac field. Zero, second, fourth har- 
monic response of CMF versus k = c/a for different intrinsic nonlinear characteristics. Parame- 
ters: //o = A*2 = 1, A*i = 2,p = 0.1, a — 0.95351. 

Fig.5. (Color online) Cases of prolate spheroidal particles: ac and dc field. Zero, fundamental, 
second, third, fourth harmonic response of CMF versus k — c/a for different intrinsic nonlinear 
characteristics (a)-(e). Parameters: ixq = iJ,2 = I, Hi = 2, p = 0.1, a = 0.1, Hdc = VS^. (f) show 
the nonharmonic response of CMF. 

Fig.6. (Color online) Cases of oblate spheroidal particles: ac and dc field. Zero, fundamental, 
second, third, fourth harmonic response of CMF versus k — c/a for different intrinsic nonlinear 
characteristics (a)-(e). Parameters: //q = A*2 = 1, A*i = 2,p = 0.1, a — 0.95351, H^c — (f) 
show the nonharmonic response of CMF. 

Fig.7. (Color online) Zero, second, fourth harmonic response of CMF versus volume fraction p 
for different intrinsic nonlinear characteristics in cases of prolate((a)-(c)), oblate((g)-(i)) spheroidal 
particles and spherical particles((d)-(f)) under ac field. Parameters: A*o = A*2 = l,A*i = 2, A; = 
1.2, a = 0.95351. 
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FIG. 1: Jian, Gao, and Huang 



16 



-0.26 ■ 
-0.28 ■ 
-0.30 ■ 
-0.32 ■ 
-0.34 ■ 
-0.36 ■ 
-0.38 ■ 
-0.40 ■ 



bet 



:bcc 



fee 



-•-^H ac=0.3 
-•--^H'ac=0.6 
■•■ ■ 5H'ac=0.9 



(a) 



0.6 0.8 



1.0 



1.2 1.4 

q 



1.6 



1.8 2.0 



35 
30 
25 
20 
15 
10 
5 

-5 
-10 
-15 



1 ' 1 


1 1 


• \ • \ • \ • <a 
















p 

• • 


m 








p 








ft 








• 








• 






(b) : 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



1.0012- 



1.0010- 



1.0008- 



1.0006- 



1 .0004 - 



1.0002- 



1 .0000 - 



(c) 



J -r ••••••••• 



0.6 



0.8 1.0 



17 



1.2 1.4 

q 



1.6 



1.8 2.0 




T ' I I r 



12- 



10- 



8- 



4^ 

2- 
0- 



— 1 — 

10 



"T — 

12 



(b) 



"T — 
14 



— 1 — 

16 



1.0030 



1.0025 



1.0015 



1.0010 



1.0005 



1.0000 



0.9995 




18 



1 ■ 1 ■ 1 


1 1 

EH^ac=0.3 - 




- - - v-\-\^ar.—n R 

— — — c.r\ ac— U.D 




£H'ac=0.9 _ 






\ / 


(prolate) , 




(a) 



1.20 

k 






0,5 1,0 1.5 2,0 2,5 3,0 3,5 4,0 4,5 5,0 






25 -| 




20- 




15- 




10- 




5- 




0- 






+ 


-5- 








-10- 




-15- 




-20- 




-25 - 













\ / ■ ; 
V ■ .' 

(f) ' 



FIG, 5: Jian, Gao, and Huang 



20 















Y 




V 


5H^ac=0.3 




!;H'ac=0.6 " 




5H'ac=0.9 




(Oblate) 
(a) 





10 12 14 16 




10 12 14 16 







■ / 






(f) 



2 4 6 8 10 12 14 16 
k 



FIG. 6: Jian, Gao, and Huang 



21 








1.0D1?- 






l.ODOS- 






1.0004- 




(b) 







0.04 0.06 




(sphere) 
(d) 

06 D.oe 



06 D.08 



0.02 0.04 0.06 





0.02 0.04 0.06 0.08 



FIG. 7: Jian, Gao, and Huang 



22 



Nonlinear response of the magnetophoresis in inverse ferrofluids 

Y. C. Jiani'2, Y. Gao^ and J. P. Huang^* 
^Surface Physics Laboratory (National Key Laboratory) and 
Department of Physics, Fudan University, Shanghai 200433 
^Department of Physics, National Tsing-Hua University, Hsinchu300 

(Dated: February 6, 2008) 

Abstract 

Taking into account the structural transition and long-range interaction (lattice effect), we resort to the 
Ewald-Kornfeld formulation and developed Maxwell-Garnett theory for uniaxially anisotropic suspensions 
to calculate the effective permeability of inverse ferrofluids. And we also consider the effect of volume 
fraction to the magnetophoretic force on the nonmagnetic spherical particles submerged in ferrofluids in 
the presence of nonuniform magnetic field. We find that the coupling of ac and dc field case can lead to 
fundamental and third harmonic response in the effective magnetophoresis and changing the aspect ratio in 
both prolate and oblate particles can alter the harmonic and nonharmonic response and cause the magne- 
tophoretic force vanish. 

PACS numbers: 83.80.Gv, 77.22.Gm, 41.20.Cv, 82.70.-y 



* Electronic address: jphuang@fudan.edu.cn 
Typeset by REVTbX 1 



I. INTRODUCTION 



Magnetic particles are used to label and manipulate biomaterials such as cells, enzymes, pro- 
teins, DNA, to transport therapeutic drugs in the applications such as bioseparation, immunoas- 
says, drug delivery, and to separate red and white blood cells yj]. Magnetizable particles expe- 
rience a force in a nonuniform magnetic field caused by the magnetic polarization. This phe- 
nomenon, called magnetophoresis (MAP) has been exploited in a variety of industrial and com- 
mercial process for separation and beneficiation of solids suspended in liquids [i2fl. Virtually, 
almost all practical magnetic separation technologies exploit the magnetophoretic forces. Such as, 
belt and drum devices, high-gradient magnetic separation (HGMS) systems [3] etc.. 

In recent years, inverse ferrofluids with nonmagnetic colloidal particles (or magnetic holes) 
suspended in ferrofluids have drawn much attention for their potential industrial and biomedical 
applications si 6[ 0]. The size of the nonmagnetic particle (such as polystyrene particle) is 
about l~100/im. Ferrofluids, also know as magnetic fluid, is a colloidal suspension of single- 
domain magnetic particles, with typical dimension around lOnm, dispersed in a liquid carrier 
Because the external magnetic field will induce magnetic dipole moment in the nonmagnetic 
particles, when the field intensity increases, the inverse ferrofluids will solidify to crystals due to 
dipole interactions, which is similar to ER and MR behavior. Since the nonmagnetic particles are 
much larger than the magnetic fluid nanoparticles, the magnetic fluid is treated as a one-component 
continuum with respect to the larger nonmagnetic particles [8,1^. Recently, Jian, Huang and Tao 
first theoretically suggested that the ground state of inverse ferrofluids is body-centered tetragonal 
(bet) latticejiol] , just as the ground states of electrorheological fluids and magnetorheological 



fluids lU 



1211 ■ Motivated by this study, we propose an alternative structural transition from the 



bet to the face-centered cubic (fee), through the application of external magnetic field. In this 
work, we will investigate the magnetophoretic force on the particles during the phase transition in 
the crystal system considering the effects mentioned above. 

Finite-frequency responses of nonlinear dielectric composite materials have attracted great at- 
tention in both research and industrial applications during the last two decades [|23|] . In particular, 
when a composite containing linear permeability particles embedded in a nonlinear permeability 
host medium is subjected to a sinusoidal alternating magnetic, the magnetic response w ill g ener- 



ally consist of frequencies of higher order harmonics under ac field 



26L 



27 



29 



sob. The 



main aim of the present paper is to study the the effects of the geometric shape and volume fraction 



2 



of the particles on the nonlinear ac responses (harmonics) of the crystal system, then the magne- 
tophoresis of the particles, considering the local lattice effect. We shall use the Edward- Komf eld 



formulation 114 , 



15i 



1811 to derive the local magnetic fields and induced dipole moments in inverse 



ferrofluids and then perform the perturbation expansion method to obtain the fundamental and 
higher harmonics in two cases: single ac and dc-ac magnetic field. 

Then, we apply a nonuniform magnetic field to investigate the magnetophoresis of the nonmag- 
netic particles submerged in the ferrofluids, by taking into account the effect of structural transition 
and long-range interaction (lattice effect). And we also consider the effect of volume fraction to 
the magnetophoretic force exerted on the nonmagnetic particles. 

This paper is organized as follows. In Sec. II, we introduce the nonlinear characteristics in 
the ferrofluids and the theory of magnetophoresis. By using the Ewald-Komfeld formulation and 
the well-known developed Maxwell-Garnett theory for uniaxially anisotropic suspensions, we cal- 
culate the effective permeability of the inverse ferrofluids for different harmonic response under 
time-varying magnetic field in Sec. III. Then we investigate the nonlinear response of magne- 
tophoretic force on the nonmagnetic particles in ferrofluids. This paper ends with a discussion and 
conclusion in Sec. IV. 
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II. THEORY AND FORMALISM 



A. Nonlinear characteristics in the ferrofluids 

In the standard literature, we often consider the magnetic induction B to be proportional with 
the field strength Hq, and their relation is written as B = /xHo, where /i is the linear permeability. 
However, the simple relation is only valid at weak intensities. When in the real situation, especially 
under a stronger field, nonlinearities are introduced as B = yuHo + ,^|//opHo + r]\Ho\^Ho + 
where ^,r] are the third-order and fifth-order nonlinear coefficient and 771-^01"^ ^ ^l-^^oP ^ A*- 
Here we assume the nonlinearity is not strong and consider only the lowest-order nonlinearity for 
simplicity. Thus after dropping the symbol of absolute value in the equation above, the nonlinear 
permeability jl is obtained: 

/i = /ii + eii^o (1) 

For inverse ferrofluids, the peameability fx of host medium ferrofluid can be derived in the form 
as /i = yUi + ^iHq when the magnetic particles in the ferrofluid is nonlinear to the external field, as 
analyzed below. As shown in Fig.(l), considering the magnetic particles with nonlinear jlp in the 
ferrofluid , the local magnetic field Hi inside the particle induced by external field Hq is give by 

Hi = Ho (2) 

Here fi j represents the host medium permeability in the ferrofluid, and the nonlinear perme- 
ability fip of the magnetic particle can be expressed as jlp = Hp + ^ Hf with ^ being the nonlinear 
coefficient. For the spherical shape particles, the effective permeability yU/ e of the whole ferrofluid 



can be determined by the Maxwell- Gamett theory [|22h 



where / denotes the volume fraction of the magnetic particles in ferrofluid. Applying Eq.(2) 
and (3) and performing the Taylor expansion by taking ^Hf as a small perturbation, we can obtain 



^^''^ (l-/K + (2 + /W ^ [{l-f)^^p + i2 + f)^,f]{^^p + 2^^fr^ 
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Hence the nonlinear permeability of ferrofluid can be expressed as 

/i = yU/,e = /ii + ^iHq (5) 

in which the constants fx, = A'/[(i+2/)/^p-2Cf-i)M/] ^ _ J,'Jn\( 4-9 ^2 + 0[e] We 

emphasize that when the magnetic particles form chains under strong external field, Eq. (3) is sup- 
posed to modify by introducing demagnetization factor |3^, but the expression of jl will remain 
unchanged by appropriately choosing fii and ,^1. 

B. Local Magnetic Field and Induced Dipole Moment 

When the inverse ferrofluids system forms crystal under external field, there are two dominant 
factors effecting its effective permeability: the local field of crystal structure and the geometrical 
shape of the particles which makes up of the crystal. Below we will first investigate the local field 
factor by performing an Ewald-Kornfeld formulation so that the structural transition and long- 
range interaction can be taken into account explicitly. Considering the ground state of inverse 
ferrofluids as a body-centered-tetragonal (bet) lattice, which can be regarded as a tetragonal lattice, 
plus a basis of two particles each of which is fixed with an induced point magnetic dipole at its 
center. One of the two particles is located at a corner and the other one at the body center of the 
tetragonal unit cell, as shown in Fig.(l). Its lattice constants are denoted by a{= h) = ^q^^ and 
c= along x{y) and z axes. In this case we identified a different structure transformation from 
the bet to fee to bee by changing the uniaxial lattice constant c under hard-sphere constraint. As q 
varies the volume of unit cell remains f c = • In this way, the degree of anisotropy of tetragonal 
lattice is measured by how q is deviated from unity and the uniaxial anisotropic axis is along the z 
axis. In particular, q = 0.87358, 1 and 23 represent the bet, bee and fee lattice, respectively. 

In case of an external magnetic field Hq along the z axis, the induced dipole moment m are 
perpendicular to the uniaxial anisotropic axis. The local field Hr at the lattice point R=0 can be 
determined by using the Ewald-Kornfeld formulation 



In this equation, 71 and 72 are two coefficients, given by, 

erfc(r/r) 2r] / 2 2n rn\ 
liyr) = o — + exp{-r] r ) (7) 
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3erftM + ,V ^ % ) ( V) (8) 

where erfc(r7r) is the complementary error function, and r/ is an adjustable parameter making the 
summations in real and reciprocal lattices converge rapidly. In Eq. (7), R and G denote the lattice 
vector and reciprocal lattice vector, respectively, 

R = C,{q~^lx + q~^my + qnz) (9) 



G = —{q^ux + q^vy + q ^wz) (10) 
where, Z, m, n, m, f , w are integers. In addition, Zj and Rj are respectively given by. 



j - 1 j - I 
Zj = I — , Rj = \R — {ax + ay + cz)\ (11) 



and n(G) = 1 + exp[i{u + v + w)/tt] is the structure factor. The local field will be computed 
by summing over all integer indices (j, /, m, n) ^ (1, 0, 0, 0) for the summation in real lattice 
and {u, V, w) ^ (0, 0, 0) for the summation in the reciprocal lattice. Here we consider an infinite 
lattice. For finite lattices, one must be careful about the effects of different boundary conditions 
I19I1 . Because of the exponential factors, we may impose an upper limit to the indices, i.e., all 
induces ranging from —L to L, where L is a positive integer. 
Now, let us define a local field factor a, 

a = (12) 

An m 

It is worth noticing that a is a function of single variable q. Also, there is a sum rule 2a± +a\\ 



=3 [2^|21I], where «|| denotes the local field factor in longitudinal field case, while a± denotes 
the local field in transverse case. Here the longitudinal (or transverse) field case corresponds to 
the fact the Hq field is parallel (or perpendicular) to the uniaxial anisotropic z axis. 



C. Nonlinear AC response in the Cases of Spheroids 

Now we are in the position to take into the account the geometrical shape effect of the particles 
by investigating nonspheres: prolate and oblate. In MR(ER or inverse ferrofluids) fluid, the shear 
stress and permeability(permittivity) will be enhanced by changing the ratio of the ellipsoidal 

6 



particles in the chain CSll such as p-phenylene-2, 6- benzobisthiazole polyparticle. This induced 
polarization can be described in demagnetizing factor, which includes the three principal axes of 
the spheroid. For a prolate spheroid(a = b < c), the longitudinal demagnetizing factor gi along 
the major axis c is given by gi 



fcln(k+v^k^) 

(fc2-i)i 



fc2-l 



with A; = ^; for oblate case(a < b = c). 



9i 



fc2 arcsin V^'^^"^ 

1 ^f^^ — ^ 1^ ■ There is a sum rule for the longitude and transverse demagnetizing 



factor [122, 



3211 that gi + 2gt = 1. For prolate spheroid, because it can easily align with its major 



axis(c) along the direction of the external field, while the dipole moments of nonmagnetic particles 
are oriented opposing the external field, here we only focus on the longitude case. It is also worth 
noticing that the spherical shape of particles is included when gi = gt = \. 



N 



1122L 



ow We can use the developed Maxwell-Gamett theory for uniaxially anisotropic suspensions 



3l|] considering both the local lattice and geometrical shape effect to derive the effective 



permeability [l32ll 



Similar to the technique applying in Eq.(4), combining Eq.(l) we admit 



gi{fio - fi) 



(13) 



HilaHi + gijxa - 3p)y] 
+ 9ixy] 



afif + gfxy^[xa — 3p] + 2agifiixy 3(?f/iQ/i^pa;^^if, 



a[y9i - i^iVii^i + gixvY 



(14) 



with X = p — l,y = fii — fiQ. 



In the presence of external oscillating time- varying magnetic field cite jlTll . the magnetic 
particles will have nonlinear characteristics and the inverse ferrofluids system will induce har- 
monic response at different frequencies without phase transition lllin . In the experiment, the 
second term or higher order terms on the right-hand of Eq. (14) can be obtained using 
mixed-frequency measurements jisl . In this paper, we will focus on the characteristics of the har- 
monics of effective permeability for the inverse ferrofluid and their nature on the magnetophoresis 
when particles put in: i^o-, Hu:-, l^2w-- In case of an external sinusoidal magnetic field(ac) of the 
form 



-o[e 



Ho = HacSinuot (15) 

Applying Eq.(15) into Eq.(14), we can obtain zero, second and fourth order harmonics for the 
effective permeability: /Xe = /^o + /X2tjCos2a;t + /i4^cos4a;t. Here we give the expression for n^^: 

7 



r2 Sgffi^px 
'8a[gixy + giy] 

On the other hand, under the application of a dc and ac field, namely. 



Hq = Hdc + Hacsmut (17) 

Using the same technique we can rewrite the effective permeability as /ig = /io + /io^sincut + 
yU2a;Cos2ci;t + yUa^sinStut + yU4^cos4co't, for extracting different harmonics. 

D. Magnetophoresis in inverse ferrofluids system 

Now we will investigate the magnetophoretic behavior when particles are put in inverse fer- 
rofluids. In a standard case, a magnetically polarizable object will be trapped in a region of a fo- 
cused magnetic field, or nonmagnetic particles in a magnetizable liquid(ferro fluid), provided there 
is sufficient magnetic response to overcome thermal energy and the magnetophoretic force ||2|]. 
Here we consider the nonmagnetically homogenous particle in the nonlinear inverse ferrofluids 
for simplicity, and the magnetic field is assumed to be dc and hysteresis is ignored. For a magneti- 
cally linear particle under magnetophoresis, the effective magnetic dipole moment vector induced 
inside takes a form very similar to that for the effective dipole moment of dielectric paricle 161. 
m = ATiRlj^jbHo, where Rejf is the effective radius of spherical or spheroidal particle and b is the 
clausius-Mossotti factor(CMF) along the direction of external field. Here we focus on the magne- 
tophoresis of spherical particle in inverse ferrofluids, and CMP b is written as: b = ^'^^^^ where 
Ii2 is the permeability of nonmagnetic particle, and jie is the permeability of the host medium: 



inverse ferrofluids. It should be noted, that in a general case, for spheroidal particle, CMF b is 

1 A'2-A'e 
3 Me+gi(M2-Me) ■ 



written as \ J^'^, ^1 — The force exerted by a nonuniform magnetic field on the dipole can be 



obtained: 

— * 

Fdipole = fJ'em ■ VHq (18) 

Combining the effective magnetic dipole moment, the magnetophoretic force exerted on the 
spherical particle in a nonuniform magnetic field can be written as, 

Fmap = 2Tt^ieR'bVHl (19) 
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Thus the magnetophoretic force on particles are proportional to the particle volume, the polar- 
ization difference and distribution of the geometric field gradient. There are two kinds of magne- 
tophoresis depended upon the relative magnitudes of ^2 and /ig. Particles are attracted to magnetic 
field intensity maxima and repelled to the field-generating electrodes when 6 > as positive mag- 
netophoresis and negative magnetophoresis corresponds to 6 < 0. Based on the harmonics of 
effective permeability, the expressions for the CMF can explicitly be obtained when the lattice 
structure of the host medium is changed, as shown in Fig. 2. The degree of anisotropy q plays 
an important role in determining the local field factor a. There is a plateau during the structure 
transformation on the relationship between q and a [|29D, accordingly we can also predict similar 
behavior in CMF, especially in Fig. 2(b). It is found that increasing q causes both zero and fourth 
order harmonic CMF decreases in Fig. 2(a) and Fig. 2(c). Note that when q changes from 0.6 to 
0.8 all CMF have great decreasing and CMF in Fig. 2(c) changes its sign when the nonlinearity is 
the highest, thus demonstrating the sensitivity of CMF when the system is under structure changes. 
Although the fourth-order harmonics of the local magnetic field and induced dipole moment seems 
to be complicated, since the strength of the nonlinear polarization of nonlinear materials can be 
reflected in the magnitude of the high-order harmonics of local magnetic field and dipole mo- 
ment, the effective nonlinear part in the effective nonlinear magnetic constant can be expected to 
determine the magnitude of the fourth-order harmonics. 



ni. NUMERICAL RESULTS 



In the preceding section we will calculate the effect of spheroidal shape and volume fraction of 
nonmagnetic particles which forms bet lattice structure in the ferro fluids, on the harmonic response 
of CMF. For a special case, when we would like to conclude the nonharmonic response situation, 
the external field Hq is constant and not vary with time namely, it can be easily obtained without 
tedious calculation by setting cut = | for ac field in Eq. (15). Thus only the zero harmonic 
contributes. Similar conclusions can be obtained that only hQ + h^ — b^^ contributes in a dc and ac 
field case. 

In Fig. 3 and Fig. 4 we show the harmonic response of CMF for different shapes in ac field: 
prolate and oblate. The zero(or nonharmonic) response of CMF exhibits strong sensitivity to the 
prolate demagnetization factor(the shape) and some peaks are obtained, while for oblate CMF bo 
have slight change compared(Fig. 4). This can be briefly explained through self-consistent ap- 
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proach 113 311 by the spectral representation separating the material parameter s = {1 — hq/ fi) from 
the geometrical parameters. Without considering the crystal lattice effect, the effective permeabil- 
ity can be derived as /ie = — g-g;('i-p) ]- It is clear that decreasing gi leads to increasing /ig, 
while for prolate(oblate) case, increasing aspect ratio k causes increasing(decreasing) of gi. In 
the meantime, because of the existence of weak nonlinearity, the nonlinear polarization has just 
a perturbation effect, which can be neglected when compared to the linear part. Furthermore, we 
can predict that one or two shape-dominant crossover ratio at which there is no net force on the 
particle in the magnetophoresis. In Fig. 3(a), it is interesting that the crossover ratio is mono- 
tonically increasing as the nonlinearity increases. Fig. 4 shows that, as nonlinear characteristics 
increase, the harmonics response of CMF are caused to increase accordingly for oblate particles. 
We emphasize to point out that, although prolate and oblate particles can both exist in different 
system state, the alignment of one particle can only be along the longest axis(prolate) when the 
particle is homogenous under the external magnetic field. 

Fig. 5 and Fig. 6(a)-(e) show different harmonic response of CMF for prolate and oblate 
particle in dc and ac field. Comparing with in single ac field, the response nature of CMF keeps 
unchanged for the same harmonics. So the introduction of coupling between dc and ac field add 
more odd harmonic response. As discussed above the nonharmonic response of CMF can be 
obtained as + h^^ — b^^j, which is shown in Fig. 5(f) and Fig. 6(f). It is also found that lower 
harmonic response of CMF can be several orders of magnitude larger than the higher one. 

Fig. 7 displays for prolate, oblate spheroidal and spherical particle, the effect of volume fraction 
of magnetic particles in the ferrofluids on the harmonic response of CMF. It is observed from 
Fig. 7(a) that the nonlinearity for nonharmonic response has little effect on CMF. The reason is 
that in the system of interest only ferrofluids are assumed to be nonlinear while the particles are 
linear. As the volume fraction of particles increases, the nonlinear component within the system 
have slight change naturally, yielding unchanged CMF accordingly. The second and fourth order 
response of CMF in Fig. 7(b) and (c) shows great variety when the volume fraction changes. This 
enhancement of the harmonics may arise from the mutual interaction between the dipole moments: 
when increasing the volume fraction, the mutual interaction becomes strong and leads to shift of 
harmonic response. It is observed that for prolate particles, the nature of zero-order response of 
CMF is greatly different from that of oblate spheroidal and spherical particles, comparing Fig. 
7(a) with Fig. 7(d) and Fig. 7(g). 
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IV. DISCUSSION AND CONCLUSION 



In this work, we have investigated both the effect of lattice structure of crystal and the geomet- 
rical shape of particles embedded in the ferrofluids on the nonlinear response of magnetophoresis 
using the perturbation method. The nonlinear response of CMF under the structure transition is 
studied in detail. It is thus possible to perform a real-time monitoring of structure transforma- 
tion. We find that the coupling of ac and dc field case can lead to fundamental and third har- 
monic response in the effective permeability constant and the CMF respectively. We also find that 
change the aspect ratio in both prolate and oblate particles can alter the harmonic and nonhar- 
monic response and thereby causing the magnetophoretic force vanish. Furthermore, by taking 
into account the local-field effect from the mutual interaction, we examine the CMF response in 
various volume fractions, and find that increasing volume fraction can result in the enhancement of 
the strengths and the difference between the different nonlinearities and characteristic harmonics. 
Thus, it is possible to obtain a good agreement between theoretical predictions and experimental 
data by suitable adjustment of both the geometric parameters(for example, the particle shape) and 
the physical parameters(for example, the temperature, component ratio in the system, saturation 
magnetization, initial magnetic susceptibility). On the contrary, such fitting is useful to obtain the 
relevant physical information of magnetophoresis in inverse ferrofluid. 

We demonstate theoretically the shape effect of the effective magnetic constant in the inverse 
ferrofluids system. It is in fact that the magnetic particles in the ferrofluid are inhomogeneous 
in the real case. A first-principle study of the dielectric dispersion of inhomogeneous suspension 



is presented 1291 13411 . and can be used to analyze the magnetic dispersion in the dilute limit, in 



which the permeability of the ferrofluid is dependent on the frequencies of external magnetic field. 



Also when the inhomogeneity is simp 



ified as the shelled spheroidal particles in the ferrofluid 



system BSll . a self-consistent method Oal can be applied for the three-component composites. 



The magnetic anisotropy as well as geometric anisotropy can both be taken into account, and will 
be reported elsewhere. 

The nonlinearity under our consideration is weak, which is common in real situations under 
moderate fields, while in case of strong nonlinearity, the perturbation approach is no longer valid 
and should be adopted using self-consistent method instead [|37|l . And in the above discussion, we 
consider the longitude case under the external field, it can also be found constricted by the relation 
gi + 2gt = 1, the harmonic response of CMF have opposite behavior. 
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Figure Captions 

Fig.l. (Color online) Schematic graph showing the inverse ferrofluid system, in which the crys- 
tal consisted of nonmagnetic particles(permeability is denoted by hq) embed in the ferrofluid(/io). 
Ferrofluids contain nanosize ferromagnetic particles(/ip) dispersed in a carrier fluid(^/). 

Fig.2. (Color online) The zero, second, fourth harmonic response of CMF versus degree of 
anisotropy q for different intrinsic nonlinear characteristics in the prolate case. The bet, bcc and 
fee lattices which are respectively related to g = 0.87358, 1.0,25 are shown(dot lines). Parameters: 
A*o = A«2 = = 2,p = 0.1, A; = 1.2. 

Fig.3. (Color online) Cases of prolate spheroidal particles: ac field. Zero, second, fourth 
harmonic response of CMF versus k — c/a for different intrinsic nonlinear characteristics. Pa- 
rameters: /Uo = = 1, yUi = 2, p = 0.1,0; = 0.95351. 

Fig.4. (Color online) Cases of oblate spheroidal particles: ac field. Zero, second, fourth har- 
monic response of CMF versus k = c/a for different intrinsic nonlinear characteristics. Parame- 
ters: //o = A*2 = 1, A*i = 2,p = 0.1, a — 0.95351. 

Fig.5. (Color online) Cases of prolate spheroidal particles: ac and dc field. Zero, fundamental, 
second, third, fourth harmonic response of CMF versus k — c/a for different intrinsic nonlinear 
characteristics (a)-(e). Parameters: ixq = iJ,2 = I, Hi = 2, p = 0.1, a = 0.1, Hdc = VS^. (f) show 
the nonharmonic response of CMF. 

Fig.6. (Color online) Cases of oblate spheroidal particles: ac and dc field. Zero, fundamental, 
second, third, fourth harmonic response of CMF versus k — c/a for different intrinsic nonlinear 
characteristics (a)-(e). Parameters: //q = A*2 = 1, A*i = 2,p = 0.1, a — 0.95351, H^c — (f) 
show the nonharmonic response of CMF. 

Fig.7. (Color online) Zero, second, fourth harmonic response of CMF versus volume fraction p 
for different intrinsic nonlinear characteristics in cases of prolate((a)-(c)), oblate((g)-(i)) spheroidal 
particles and spherical particles((d)-(f)) under ac field. Parameters: A*o = A*2 = l,A*i = 2, A; = 
1.2, a = 0.95351. 
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FIG. 1: Jian, Gao, and Huang 
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